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Poaching is a widespread and well-appreciated problem for the conservation of many threatened species. 
Because poaching is illegal, there is strong incentive for poachers to conceal their activities, and con- 
sequently, little data on the effects of poaching on population dynamics are available. Quantifying 
poaching mortality should be a required knowledge when developing conservation plans for endangered 
species but is hampered by methodological challenges. We show that rigorous estimates of the effects of 
poaching relative to other sources of mortality can be obtained with a hierarchical state -space model 
combined with multiple sources of data. Using the Scandinavian wolf (Canis lupus) population as an 
illustrative example, we show that poaching accounted for approximately half of total mortality and 
more than two-thirds of total poaching remained undetected by conventional methods, a source of mor- 
tality we term as 'cryptic poaching'. Our simulations suggest that without poaching during the past 
decade, the population would have been almost four times as large in 2009. Such a severe impact of 
poaching on population recovery may be widespread among large carnivores. We believe that conser- 
vation strategies for large carnivores considering only observed data may not be adequate and should 
be revised by including and quantifying cryptic poaching. 
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1. INTRODUCTION 

The illegal killing of animals, hereafter poaching, threa- 
tens the viability of many species worldwide [1-5]. 
Because of their characteristic low densities combined 
with their slow rates of population growth, top predators 
are particularly vulnerable to effects of poaching. Almost 
all large carnivore species have endured a long history of 
human persecution and have been eradicated from sub- 
stantial parts of their historical ranges [6]. Although 
most species of large carnivores are now legally protected, 
poaching remains a widespread problem for their conser- 
vation [6]. Some species are commercially poached for 
pelts or body parts used in traditional medicine [7], but 
many are killed because of conflicts with human interests, 
such as competition for game, depredation of livestock 
and threats to human safety [8]. It follows that dealing 
with poaching mortality often emerges as a required con- 
dition for the restoration, conservation and sustainable 
management of large carnivore populations. 
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A near universal problem with understanding poaching 
is the absence of rigorous estimates of its effects relative to 
other sources of mortality [1]. There are several recent 
attempts to assess the extent, mechanisms and effects of 
poaching [2,9-12] but remarkably little quantitative 
data exist, although new methods to measure its extent 
have recently been developed [13]. One obvious reason 
for the absence of data is methodological. The most 
reliable method of quantifying causes of mortality in 
populations of large wild mammals is to observe their 
fates over time using radio-tracking [14]. However, 
when a radio-collared animal is poached, there is a high 
probability that the poacher promptly destroys the trans- 
mitter and hides (or consumes) the carcass, leaving the 
researcher with a lost radio contact without known 
cause [15]. Treating cases of lost radio contact in a sur- 
vival analysis based on radio-tracking is not a trivial 
problem, especially not for such 'poaching-prone' animals 
as large carnivores. One can never exclude the possibility 
that a certain proportion of animals with lost radio con- 
tact in fact died from poaching that cannot be verified. 
We define this unobserved source of mortality as 'cryptic 
poaching'. Estimating a quantity in ecological processes 
that is not amenable to direct observation is feasible 
with hierarchical models because these models allow 
multiple sources of data to inform estimates of model 
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parameters, including unobservable ones [16]. These 
sources of data can include observations on state variables 
from monitoring studies as well as direct estimates of 
observable parameters from detailed studies of processes. 
In this paper, we used a decade (1999-2009) of popu- 
lation census, radio-tracking and recruitment data of the 
Scandinavian wolf population combined with a Bayesian 
state -space hierarchical population model to show that 
poaching has drastically slowed down the recovery of 
this population. 



2. MATERIAL AND METHODS 

(a) General approach 

Between December 1998 and April 2009, we radio-marked 
104 wolves in Scandinavia constituting between 10 and 15 
per cent of the population, among which we had 26 verified 
mortalities. We used radio-tracking data to compute three 
cause-specific mortality rates based on 21 cases of non- 
poaching (seven natural deaths such as age and disease, 
five traffic mortalities and nine cases of legal control), five 
cases of verified poaching, and finally, 18 cases of cryptic 
poaching (not included in the 26 verified mortalities). We 
considered a wolf as having been cryptically poached or ver- 
ifiably poached according to criteria explained below. 
However, we could not obtain a robust estimate of cryptic 
poaching because we never found the supposedly dead 
wolves. We circumvented this obstacle by fitting a hierar- 
chical state -space model to another dataset, a decade-long 
time series of population size and number of reproductions. 
In particular, we investigated whether the non-poaching 
mortality and the verified poaching rates would be large 
enough altogether to explain the observed population 
trends, or, on the contrary, if an additional source of 
mortality was needed to fit the longitudinal data the best. 

(b) Criteria for cryptic and verified poaching 

Cryptic poaching was defined based on four criteria (with 
either all of criteria 1-3 or criterion 4 alone satisfied): 

1 . Sudden loss of radio contact with no indication of trans- 
mitter failure (more than half of the expected battery 
life-time remaining). 

2. At least two aerial searches over a much larger area than 
the wolf territory were performed without further contact 
with the collared individual. 

3. The individual was resident and repeated snow-tracking 
within the territory, in combination with the collection 
of scats and subsequent DNA analyses of multiple 
faeces confirmed that this individual was no longer 
present within the pack territory. 

4. Radio contact was lost and special circumstances strongly 
indicated that poaching was the most plausible expla- 
nation. This applied only for two cases where police 
reports confirmed that people had attempted to poach 
wolves. 

Wolves not satisfying these criteria were censored at the date 
of lost contact. 

Verified poaching was defined based on two criteria 
(enough if one criterion is satisfied): 

1 . The body was recovered and the necropsy showed that a 
human deliberately killed it outside a legal hunt. 



2. Wolf tissue (skin or muscle) determined by DNA analysis 
to originate from one of the radio-collared wolves was 
found in possession of a person that could not explain 
how he had acquired it and was later convicted at a 
court for this illegal possession. 

(c) Hierarchical model 

To estimate the posterior distribution of the true size of the 
population, we composed process and observation equations. 
The process equation was 

IX, = log[iV,_i(l -m-v-c) + IRt-i] 
and 

Nt - lognormal(^,, aproc), 

where is the deterministic prediction of the median wolf 
population size at time r, is the true population size at 
time t, cTpj-oc is the standard deviation of the true population 
size on the log scale, m is the mortality rate from all causes 
except poaching, v is the verified poaching rate, c is the cryp- 
tic poaching rate, / is the per pack recruitment rate and is 
the number of reproductions at time t. The process equation 
was linked to data using the observation equation 



Xt - gamma (t^t, A) 

and NobSt ~ Poisson(At), 

where ATobs^ is the observed population size at time r, cr|^obs is 
the estimate of the error of observation of the population size. 
This formulation views the count data hierarchically — the 
mean observed count of wolves at time t is Poisson distribu- 
ted with mean and this mean is drawn from a gamma 
distribution with mean equal to the prediction of the process 
model and a standard deviation for observation error. We 
chose this approach because it allows the uncertainty in the 
data model to be larger than the variance of the Poisson par- 
ameter Af. The approach is the same as assuming that the 
count data follow a negative binomial distribution, but offer 
computational advantages [17]. We did not include density- 
dependence in our hierarchical model because there is plenty 
of space and wild ungulate prey for larger wolf populations 
on the Scandinavian Peninsula, and both Sweden and 
Norway have some of the highest moose-wolf ratios in the 
world [18]. 

(d) Data and model priors 

Estimates of total population size (ATobs) and number of 
reproductions (R) were obtained annually from 1999 to 
2009, using a combination of snow-tracking, radio-tracking 
and DNA analysis of scats (see electronic supplementary 
material, time-series data). 

Monitoring of pack reproductions provided informative 
prior on /. Number of pups at the age of six months was esti- 
mated from recurrent sessions of snow-tracking within 
territories (3.788 + 1.466). Shape parameters of informative 
gamma-distributed prior for litter size / were then calculated 
using moment matching (table 1) [19]. 

Using radio-tracking data, we calculated cause-specific 
mortality rates: non-poaching w = 0.148 + 0.028, verified 
poaching = 0.049 + 0.017, and suspected poaching 
c = 0.085 + 0.023, and accounting for competing risks [14]. 
Shape parameters of informative beta-distributed priors 



Proc. R. Soc. B (2012) 



912 O. Liberg et al. Quantifying cryptic poaching 



Table 1. Parameter values from radio-tracking dataset and 
prior distributions. Priors shape parameters were derived by 
moment matching of mean and s.d. values from the radio- 
tracking dataset. Note that o-proc is on a log scale. 



mean value 
from radio- 
parameter tracking s.d. prior 



m 


0.148 


0.028 


beta (23.44, 135.21)^ 


V 


0.050 


0.017 


beta (7.84, 150.56)^ 


c 


0.085 


0.023 


unif (0, 1) 


I 


3.788 


1.466 


gamma (6.67, 1.76)^ 








unif (0, 50) 


^proc 






unif (0, 25) 



^Denotes informative priors. 



for rates m and v were calculated using moment matching 
(table 1) [19]. Although radio-tracking data could provide 
information also on cryptic poaching rate c, we specified an 
uninformative prior on c to investigate whether the observed 
population trends could be well explained with no cryptic 
poaching at all. While literature refers to cause-specific mor- 
tality rates, c, m and v are in fact probabilities to die from a 
specific cause, and therefore, the most uninformative prior 
we could give to c was a uniform distribution in (0,1). 

Shape parameters for uninformative priors for cr^obs and 
cTproc were chosen subjectively to assure that all possible 
values in the posterior distribution had equal densities in 
the prior distribution (table 1). Sensitivity of posterior dis- 
tributions to the priors was tested to assure the choices of 
shape parameters were uninformative. 

(e) Monte Carlo Markov Chain inference 

We estimated the posterior distribution of each parameter by 
running Monte Carlo Markov Chains, implemented in JAGS 
[20] with R [21]. Six chains were initialized with different 
sets of parameter values chosen within biologically plausible 
bounds. After an initial burn-in period of 100 000 iterations, 
we obtained 1 000 000 iterations of each of the chains, thinning 
each by 10. We successfully checked for convergence using the 
Heidelberger & Welch [22] stationarity and half- width tests 
with the CODA package [23]. We evaluated the overlap 
between prior p{6) and posterior 7r{6\y) distributions by com- 
puting the quantity t= {\ / p{6)) I mm{ p{6) , 7r{0\y))dO [24]. 

To estimate the impact of poaching on our study popu- 
lation, we simulated population trajectories from 1999 to 
2009 using posterior distributions of parameters without 
cryptic poaching and without any poaching at all. We also 
wanted to investigate if we could differentiate an absence of 
cryptic poaching from small rates of cryptic poaching. For 
this, we considered that the simulated population without 
cryptic poaching would become our longitudinal data and 
we used the same approach to estimate the posterior distri- 
bution of cryptic poaching. This amounts to fitting a model 
to a dataset that we know has been generated with no cryptic 
poaching at all. 

3. RESULTS 

Poaching accounted for half of total mortality (5 1 %) and 
more than two-third (69%) of total poaching was cryptic. 
The median estimates of posterior non-poaching 



(0.142 + 0.027) and verified poaching (0.046 + 0.016) 
mortality rates from the model were very similar to the 
rates based on radio-tracking data (respectively, 0.148 + 
0.028 and 0.049 + 0.017; figure 1). The median estimate 
of the posterior cryptic poaching rate from the model 
(0.103 + 0.106) was also remarkably close to the indepen- 
dent estimate based on radio-tracking data (0.085 + 
0.023) but was accompanied by higher variance (figure 1). 
However, despite this variance, the data improved the 
estimate of cryptic poaching over prior knowledge. 
Overlap was large for non-poaching (r = 94%) and verified 
poaching {t= 96%) mortality rates. On the contrary, for 
cryptic poaching overlap was smaller {t= 37%) indicating 
that modelling did add information on the estimate of 
cryptic poaching. 

Our study population increased from 74 individuals in 
winter 1998/1999 to 263 in 2008/2009 (figure 2). The 
mean annual growth rate during this period was 13.5 
per cent. Assuming no (verified and cryptic) poaching 
and no density-dependence, this trajectory would have 
resulted in a median population size of 990 wolves in 
2009, i.e. almost four times larger than the one observed. 
For a population without cryptic poaching and with veri- 
fied poaching only, the trajectory would have resulted in a 
median population size of 676 wolves in 2009. When we 
considered the simulated population without cryptic 
poaching as data and used the same modelling approach 
to quantify cryptic poaching, we obtained a rate of cryptic 
poaching very close to zero (c = 0.023 + 0.03). 

4. DISCUSSION 

Here we have demonstrated a high incidence of poaching 
in a threatened wolf population, but because a major part 
of this poaching was unobserved (cryptic) and inferred 
from indirect data, its estimate is open to criticism. 
Although we cannot unequivocally prove that the 
posterior of rate c strictly includes only poaching, we 
can conclude that an additional source of mortality is 
required to explain our combined data. Because we 
could not identify any other cause of death than poaching 
that also would have resulted in a sudden loss of radio 
contact combined with no later verification through the 
continuous DNA-monitoring of the population, we 
believe that most, if not all, mortality included in this 
rate is indeed cryptic poaching. The close correspondence 
between the estimates we obtained for cryptic poaching 
rate from the model (0.103) and the independent one 
based on radio-tracking data (0.085) furthermore support 
that cryptic poaching indeed was an important mortality 
cause in our study population. The larger model estimate 
might be explained by the fact that the estimate based on 
radio-tracking data could be an underestimate. Because of 
our strict criteria for poaching, we did not classify any of 
the missing radio-collared non-resident dispersers as 
cryptic poaching, although this may have occurred in 
some cases. Our estimate of cryptic poaching received 
further support by the large gap between the simulated 
population trend without cryptic poaching and the 
observed dataset. This should convincingly reveal that 
observable mortality rates m and v cannot alone explain 
the observed population trends. 

Our results may have been severely biased if we had 
underestimated population size, because, in that case, 
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Figure 1 . Posterior (soHd black Hne) and prior (dotted black line) densities for (a) cryptic poaching rate c (posterior median = 
0.103 + 0.1 06^ shown by vertical thin line)^ (6) non-poaching mortality rate m (posterior median = 0.142 + 0.027) and 
(c) verified poaching rate v (posterior median = 0.046 ± 0.016). Overlap between prior and posterior densities is shown 
by the grey area. Parameters m and v were given informative priors based on radio-tracking data. Their posterior median esti- 
mates were very similar to rates from radio-tracking data (non-poaching mortality rate = 0.148 + 0.028, verified poaching 
rate = 0.049 ± 0.017). The prior for cryptic poaching rate was on the contrary left uninformative. Still, its posterior median 
estimate was remarkably similar to the independent estimate of cryptic poaching rate from radio-tracking data (0.085 ± 
0.023). The posterior density of cryptic poaching poorly overlapped with its prior and reveals that an unobserved source of 
mortality was present in the population. 



the parameter we attributed to cryptic poaching would 
actually have been a correction factor for our underesti- 
mate. However, the intense and continuous fine scale 
monitoring of wolves allowed us to rule out a systematic 
underestimate of population size (see electronic sup- 
plementary material, robustness of census data). Our 
results may have been equally biased if our criteria for 
cryptic poaching were inadequate. However, none of the 
18 animals classified as cryptic poaching were ever 
detected after the loss of radio contact by any of the 
survey methods used in this study (see electronic sup- 
plementary material, robustness of poaching criteria). 
Quantifying a cause-specific mortality rate based on 
unknown fates requires also excluding the possibilities 
that an animal would have remained undetected by dis- 
persing from the population. The breeding wolf 
population on the Scandinavian peninsula (Norway and 
Sweden) appears to be functionally isolated from the 
Finnish -Russian population with very little immigration 
and only one confirmed emigration recorded during the 
past decade (see electronic supplementary material, 
population isolation). 

Poaching has had a significant impact on the popu- 
lation recovery. An average annual growth rate of 
13.5 per cent is well below the typical rate of colonizing 
or recovering wolf populations [25]. Without any 



poaching, the median annual growth rate of our study 
population would have been 29.5 per cent during the 
period, i.e. more than double the observed rate and 
compatible with the fastest recovering wolf populations 
on record [25,26]. Considering that neither suitable 
habitat nor prey base are limiting factors, the population 
size in 2009 would probably have been three to four 
times the one observed (figure 2). Although the popu- 
lation has continued to grow, the decelerated growth 
rate caused by poaching is having other negative conse- 
quences. It has postponed the time when managing 
authorities can be more fiexible with permits to kill pro- 
blem individuals, causing unnecessary confiict with local 
people. Still more serious, it has aggravated an already 
bad genetic situation. The Scandinavian wolf population 
is small, isolated and facing serious genetic problems 
[27,28], and any delay in growth will accelerate inbreed- 
ing and loss of genetic variation [29] . 

Few studies of large carnivore survival based on radio- 
tracking have clearly described how they have treated cases 
of lost radio contact and made efforts to differentiate 
between possible fates of these animals. In a newly protected 
wolf population in north central Minnesota, a substantial 
proportion of lost radio contacts was assumed to be caused 
by illegal killing, and estimated to make up 70 per cent of 
total mortality rate [30]. In three different Scandinavian 
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Figure 2. Model and census estimates of the wolf population 
in Scandinavia during 1999-2009. Filled black triangles are 
census data. Squares are median of posterior distribution of 
the fitted models with its 95% credible interval shown by 
dashed lines. Circles are the median posterior distribution 
of the simulated population without poaching assuming no 
density-dependence. This reveals a decade of poaching 
scaled down population size from 990 to 263 wolves in 
2009. Lozenges are the median posterior distribution of the 
simulated population without cryptic poaching assuming 
no density-dependence. The grey area indicates the 
number of wolves theoretically lost due to cryptic poaching. 



lynx {Lynx lynx) populations, formal criteria resembling the 
ones we have used in this study were set up for validation of 
each case of lost radio contact, resulting in poaching rates 
between 32 and 74 per cent of total mortality [31]. A 
study on wolverines {Gulo gulo) also differentiated among 
lost radio contacts, and estimated that poaching made up 
60 per cent of total mortality [32], while in a study of 
Amur tigers (Panthera tigris), ratio of poaching to total mor- 
tality was 75 per cent [15] . A shared result in these studies 
was that a substantial part of the estimated poaching rates 
was made up of cryptic poaching (44-7 1 %) . Cryptic poach- 
ing was estimated to be 69 per cent of the total poaching rate 
in our study, which falls within the range of these earlier 
studies. Had all cases of lost radio contact just been 
censored from further calculations, both poaching and 
total mortality rates would have been seriously underesti- 
mated. However, as the estimates of the cryptic part of 
poaching in all these case studies were based on assump- 
tions, the degree of uncertainty in the estimates was 
unknown. By using a hierarchical model, we could combine 
multiple sources of data in a statistically coherent way and 
unobserved quantities could be estimated because of their 
interdependence with the quantities that are observed. 

We believe that the results presented above, motivate 
careful reconsideration of the extent of cryptic poaching 
in all studies of large carnivores. A recent example is the 



extensive study of survival in the newly recovered 
Northern Rocky Mountain wolf population in northwes- 
tern United States [33], where a minimum of 87 (24%) of 
the 363 dead animals was confirmed illegally killed. A 
further number of 150 animals were censored at the 
date of lost radio contact. Although the authors gave sev- 
eral arguments why it was less likely that these animals 
might have been poached, we caution that this indeed 
might be the case for a substantial part of them, especially 
considering that the study was performed in an extremely 
wolf hostile human environment. 

We have shown that the failure to include the effects of 
cryptic poaching can cause serious errors in the estimation 
of the potential rate of population growth. Because a sub- 
stantial part of poaching is often unobserved, poaching 
may be an even larger problem in wildlife conservation 
than has hitherto been assumed owing to the difficulty of 
measuring it properly. Quantifying cryptic poaching and 
its impact illustrates a challenging problem that is not 
unusual in ecology and conservation biology — the esti- 
mation of unobservable parameters with small values but 
high variance [34]. As we have illustrated here, such 
problem can be successfully addressed by combining 
multiple data in a hierarchical framework to obtain robust 
inferences. The increasing possibilities to mark many 
more individual animals at a much larger range of taxa, 
body sizes and length of tracking time [35,36] should 
make collecting individual data more feasible, and there- 
fore, our approach more widely applicable in the future. 
Our study should further reinforce the need to bring uncer- 
tainty to the centre stage of conservation studies [37] and 
illustrate how considering uncertainty affects the ability 
to manage populations [3,38]. 

All research presented complies with the current laws of the 
country in which the experiments were performed. 
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H. Ronning, R. Skyrud, T. H. Stromseth and O.-K. Steinset 
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Swedish University of Agricultural Sciences, Norwegian 
Directorate for Nature, Management, Norwegian Research 
Council, Norwegian Institute for Nature Research, Hedmark 
University College, County Governors of Hedmark and 
Varmland, Borregaard Skoger, Glommen Skogeierforening, 
Norskog, Norges Skogeierforbund, OUe and Signhild 
Engkvists Stiftelser, Carl Tryggers Stiftelse, Marie Claire 
Cronstedts Stiftelse, Swedish Carnivore Association, and 
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